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ABSTRACT 

We present the results of our optical monitoring campaign of the X-ray source 
H 0507+164, a low luminosity Seyfert 1.5 galaxy at a redshift, z — 0.018. Spectroscopic 
observations were carried out during 22 nights in 2007, from the 21 of November to 
the 26 of December. Photometric observations in the R-band for 13 nights were also 
obtained during the same period. The continuum and broad line fluxes of the galaxy 
were found to vary during our monitoring period. The R-band differential light curve 
with respect to a companion star also shows a similar variability. Using cross correla- 
tion analysis, we estimated a time delay of Teen = S.Oll^ g4 days (in the rest frame), 
of the response of the broad line fluxes to the variations in the optical continuum 
at 5100 A. Using this time delay and the width of the line, we estimated the radius 
for the Broad Line Region (BLR) of 2.53^°'g5 x 10"'^ parsec, and a black hole mass 
of 9.62l° ji^ X lO^M©. 
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1 INTRODUCTION 

Accretion of gas onto a Super Massive Black Hole (SMBH) 
in the nucleus of galaxies is believed to be the source of 
activity in Quasars and Seyfert galaxies (commonly known 
as Active Galactic Nuclei (AGNs); cf. Rees 1984). Several 
studies have suggested that the mass of the SMBH in these 
objects is correlated with the luminosity, mass and veloc- 
ity dispersion of the stellar spheroid of the galaxies (Kor- 
mendy & Richstone 1995; Magorrian et al. 1998; Ferrarese 
& Merritt 2000; Gebhardt et al. 2000; Marconi & Hunt 2003; 
Haring & Rix 2004). Such correlations may imply an evolu- 
tionary relation between the growth of the SMBH and the 
host galaxy itself (e.g. Somerville et al. 2008; Shankar et al. 
2009; Hopkins & Hernquist 2009). In order to study the de- 
pendence of the various observed phenomena of AGNs on 
the black hole mass and the cosmic evolution of the black 
holes, independent and reliable estimates of the mass of the 
black holes are required (e.g., Goulding et al. 2010; Rafter, 
Crenshaw & Wiita 2009). 

One independent method to estimate the mass of the 
black hole is using the reverberation mapping technique 
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(Blandford & McKee 1982; Peterson 1993). In the opti- 
cal bands, the continuum flux of some AGNs, is known to 
vary on timescales as short as hours (e.g., Miller, Carini & 
Goodrich 1989; Stalin et al. 2004). If the main source of 
ionization of the Broad Line Region (BLR) is the contin- 
uum itself, any variation of the continuum emission can also 
be seen in the broad emission lines. However, the variations 
in the broad line flux will have a time lag (r) relative to 
the continuum variations, which can be interpreted as the 
light travel time across the BLR. As a first approximation, 
therefore, the size of the BLR is Rblr ^ cr, where c is the 
velocity of light. Once the Rblr is obtained, the mass of 
the black hole can also be estimated, using the velocity dis- 
persion of the broad component of the emission lines, aune, 
and assuming virial equilibrium (Peterson et al. 2004; P04); 
see Peterson 2010, for a recent review). 

The reverberation mapping technique has been used to 
make estimates of SMBH masses over a large range of red- 
shift. However, because the technique is observationally tax- 
ing, as it demands an enormous amount of telescope time, 
to date the BLR radius of only about three dozen AGNs 
(Seyfert 1 galaxies and Quasars) have been determined (P04; 
Kaspi et al. 2007; Bentz et al. 2009a; Denney et al. 2009, 
2010). Nevertheless, using these estimates a correlation was 
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found between Rblr and the optical continuum luminosity 
at 5100 A (Kaspi et al. 2000; Kaspi et al. 2007; P04; Den- 
ney et al. 2009; Bentz et al. 2009b). The Rblr-\1^^_^^^^ 
relation can be considered well constrained between the lu- 
minosities lO'*'^ ere sec~^ < AL ° < 10'*^ ere sec~^. On 

5100A 

the other hand, for luminosities below lO**^ erg sec~^, only 
a handful of sources are observed, and the estimated values 
of Rblr could also indicate a flattening of the relation (see 
Fig. 2 of Kaspi et al. 2005). This flattening would suggest 
a lower limit in the possible masses of SMBHs in galaxies. 
Although recent revisions of a few sources made by Bentz 
et al. (2006) and Denney et al. (2009;2010) are consistent 
with a continuation of the Rblr. — AL_^ ° relation to lower 
luminosities, and consequently with no lower limit in the 
mass for the SMBH, the correlation is still sparsely sam- 
pled. Moreover, the Rblr — -^L^^ relation is very useful 
for estimating the SMBH masses ?rom single-epoch spectra 
and calibrating other surrogate relations used for black hole 
mass estimates (Vestergaard 2004; Shen et al. 2008). There- 
fore, estimates of Rblr for a larger number of sources are 
required. 

The extrapolation of the known Rblr ~ AL ° rela- 
tion to low lummosities suggests that the time lag between 
the variations of the broad line and that of the continuum 
will be of the order of hours to days, as compared to sev- 
eral months for high luminosity sources. Thus, monitoring 
programs of short durations, but fast sampling, are required 
to estimate the reverberation time lags for low luminosity 
sources. 

In this paper, we present the optical spectroscopic and 
photometric observations of a new low luminosity AGN, the 
X-ray source and Seyfert 1.5 galaxy H 0507-1-164. Based on 
a reverberation mapping campaign that lasted for about a 
month, during November-December 2007, we have obtained 
Rblr and estimated the mass of the SMBH. In Section 2, 
the observations and data reductions are described. The re- 
sults of the analysis are given in Section 3, and the conclu- 
sions are presented in Section 4. 



Spectrograph and Camera (HFOSCfl In imaging mode, 
only the central 2048 x 2048 pixels region of the CCD is 
used. The camera has a plate scale of 0.296 arcsecond/pixel, 
which yields a field of view of 10 x 10 square arcmin. 

2.1 Spectroscopy 

Medium resolution spectra of the nucleus were obtained us- 
ing a 11 arcmin x 1.92 arcsec wide slit and a grism. The 
spectra have a spectral range of 3800-6700 A with a reso- 
lution of ~8 A. The exposure time varied between 900 and 
1000 seconds. The spectra were reduced using standard pro- 
cedures in After bias subtraction and flat fielding, 
one dimensional spectra were extracted and calibrated, in 
wavelength using an FeAr lamp, and in fiux using various ob- 
servations of the spectrophotometric standard star Feige 34. 
Since the observed spectra were of low S/N, for further anal- 
ysis all the spectra were smoothed to a resolution of ^^15 A. 

The standard technique of spectral fiux calibration is 
not sufficiently precise to study the variability of AGNs. 
Since even under good photometric conditions the accuracy 
of spectrophotometry is not better than 10% (Shapovalova 
et al. 2008), we used a relative calibration procedure. A first 
order fiux calibration was first obtained in the normal way 
using the standard star. Then all the spectra were inter- 
calibrated relative to the spectra of one night (we choose 
the 21 of November), assuming the fiux of the narrow line 
[O in]A5007 A is constant. This is justified, because the Nar- 
row Line Region (NLR) is much more extended (of the order 
of a few hundred parsecs) than the BLR (much less than a 
parsec) and fiux variation cannot be observed in this region 
over short time scales (cf. Osterbrock 1989). Each spectra 
were scaled relative to the reference spectra using the scal- 
ing algorithm devised by van Groningen & Wanders (1992). 
This algorithm uses a chi-square, to minimize the residuals 
of the [O lll]A5007 A line after subtraction from the refer- 
ence spectrum. The mean spectrum, averaging the 22 nights 
of observations, is shown in Fig. \T\ 



2 OBSERVATIONS AND REDUCTIONS 

Using the Veron-Cetty & Veron catalogue of Quasars and 
Active Galactic Nuclei (12th Ed.; Veron-Cetty & Veron 
2006), we have compiled a list of nearby Seyfert 1 galax- 
ies, which, based on the available spectra, have a luminosity 
at A5100 A of the order of 10^^ erg sec~^ or lower. Very 
few candidates were found (mostly because of the absence 
of available spectra). The source, H 0507-1-164, that we se- 
lected for our campaign is identified in the catalogue of 
Veron-Cetty & Veron as an X-ray source, with coordinates 
02000 = 05''10'"45.5^(52ooo = 16'*29'"56", and is classified 
as a Seyfert 1.5 galaxy at a redshift of z = 0.018. 

Optical spectroscopic and photometric observations of 
H 0507-1-164 were carried out in 2007 between 21 of Novem- 
ber and 26 of December at the 2m Himalayan Chandra 
Telescope (HCT), operated by the Indian Institute of As- 
trophysics, Bangalore. The telescope is equipped with a 
2048 X 4096 CCD, coupled to the Himalayan Faint Object 



2.2 Photometry 

In parallel to the spectroscopic observations, R-band im- 
ages were also obtained with the HFOSC. Unfortunately, 
some observations turned out to be contaminated by the 
light of the extremely bright stars located near the source. 
As a consequence, our R-band photometry covers only 13 of 
the 22 nights of the campaign. The images were bias sub- 
tracted and flat fielded using IRAF packages. For the rest of 
the reduction, calibration and analysis, packages in MIDA^ 
were used. After removing the cosmic rays, profile fitting 
photometry was done using the DAOPHOT and ALLSTAR 
packages. The observed R-band frame is shown in Fig. [21 

^ http: //www.iiap.res.in/iaojifosc 

^ IRAF stands for Image Reduction and Analysis Facility and 
is distributed by the National Optical Astronomy Observatories, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc. under contract to the National Science Foun- 
dation 

Munich Image Data Analysis System; trademark of the Euro- 
pean Southern Observatory 
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Figure 1. Average spectrum of H 0507+164 obtained by averaging the observations of the 22 nights. The flux is in the unit of 
10~^® erg cm~^ s~^ A~^. The classification as Seyfert 1.5 is confirmed, as one can see both a broad and narrow component for the H/3 
hue at A4861 A(cf. Osterbrock 1989). 
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Table 1. Continuum and iT^ fluxes for the source H 0507+164 



Juhan Date 



Fx X 10-^" 

5100 A 4861 A 

-2454420 (erg/s/cm2/A) (erg/s/cm^/A) 



Am 



(mag) 



Figure 2. The R-band image of H 0507+164. The galaxy (in 
the center) and comparison star (on the top) used for differential 
photometry are indicated by circles. The field of view has a di- 
mension of 10 X 10 arcmin. North is up and east is to the left. 
Note the presence of an extremely bright star near H 0507+164. 



3 ANALYSIS 
3.1 Lightcurves 

The lightcurves of the Hp flux and the continuum at 5100 A 
were obtained using the final inter-calibrated spectra. The 
continuum flux in the rest-frame of the galaxy at 5100 A 
was obtained using the mean flux within the observed band 
from 5172 to 5200 A. 

The H/3 emission line fluxes were obtained by integral- 
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ing the emission proflle in the band spanning 4884—5012 
A, after subtracting a continuum. An average of the mean 
fluxes in the regions on the blue (4808-4852 A) and red 
(5012—5024 A) sides of the H/5 line was used as the con- 
tinuum below the line. Although the measured line fluxes 
include both the narrow and broad components, any vari- 
ation observed in the line fluxes can be attributed to the 
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Figure 3. The R-band differential lightcurve (top panel), fol- 
lowed by the lightcurves for [O in]A5007 A, H^g and continuum 
at 5100 A are plotted. The fluxes are in units of 10"^'^ 
cm-2 A-i for the [O lll]A5007 A Une, lO"!* "'^ 
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for the H/3 line and V)-^^ 
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broad component only, since the narrow component is not 
expected to vary during the period of our observations. The 
lightcurves for the continuum at 5100 A, for the H/j and 
for the [O in]A5007 A lines are shown in Fig. [S] The cor- 
responding fluxes for the continuum and H/j are listed in 
Table [T] As expected from the inter-calibration procedure, 
the lightcurve for [O in]A5007 A is nearly constant. On the 
other hand, both the continuum at 5100 A and the flux 
are observed to vary. 

The observed R-band differential instrumental magni- 
tudes between the galaxy and the comparison star (marked 
in Fig. |3]) are also given in Table [1] The R-band differential 
lightcurve plotted in Fig. [3] (top panel) closely follows the 
lightcurves of the continuum at 5100 A and H/j. The average 
flux at 5100 A is 9.03± 1.47 x 10"^® erg cm-^s"^A-\ which 
corresponds to AL •> of 3.4 ± 0.55 x 10"*^ erg s~^, for the 

cosmological parameters Ho — 70 km s~^ Mpc~^, il^ 
and Q.\ = 0.7. 



The variability of the light curves are characterised by 
the parameters, excess variance, F„ar and the ratio between 
the maximum and minimum flux of the hght curves, Rmaa: 
(Rodriguez-Pascual et al. 1997; Edelson et al. 2002). The 
continuum and H;3 line have Y^ar of 0.16 and 0.18 and Rmaa: 
of 1.9±0.020 and 1.8±0.016 respectively. These values are 
within the range of values found by other variability studies 
of AGNs (cf. P04; Denney et al. 2010). 



3.2 Time Lag 

The time lag between the variations of the continuum flux 
and the variations of the emission can be determined 
by cross-correlating the two hght curves. For cross correla- 
tion analysis, the method of interpolated cross-correlation 
function (ICCF; Gaskell & Sparke 1986; Gaskell & Peter- 
son 1987) and the method of discrete correlation function 
(DCF; Edelson & Korlik 1998) were used. Ahhough both 
ICCF and DCF methods produce similar results (White & 
Peterson 1994), the interpolation of the light curve during 
the period of gaps required by the ICCF method might not 
be a reasonable approximation to the behavior of the light 
curves. Thus the DCF method is preferable for data with 
large gaps (Denney et al. 2009). For comparison, we obtain 
the results using both the methods. 

The results of the cross correlation analysis are shown in 
Fig. [4] The cross correlation function (CCF) obtained using 
the ICCF method is plotted as a thick solid line. The auto 
correlation functions (ACFs) of the continuum at 5100 A 
and the Hp line are also shown in Fig. [4] as dashed and dot- 
ted lines respectively. For comparison the cross correlation 
function obtained using the DCF method is also plotted as a 
thin solid line. As expected, the auto correlations have zero 
time lags. On the other hand, the time lag in the cross cor- 
relation curve is clearly noticeable as an overall shift to the 
right. The position of the maximum in the CCF provides an 
estimate of the time lag between the continuum and the H/3 
line. The maximum was however determined using the cen- 
troid, which gives a better estimate for noisy CCFs, rather 
than the peak, using the formula: 



(1) 



0.3 



The estimate of the centroid includes all the points that 
are within 50% of the peak value of the CCF. Based on 
this cross correlation analysis, a statistically significant cen- 
troid and the associated uncertainty were obtained using a 
bootstrap technique that introduces effects of randomness in 
fluxes and sampling of the light curve (cf. P04). A method 
to carry out Monte-Carlo (MC) simulation using the com- 
bined effects of flux randomization (FR) and the random 
subset selection (RSS) procedures is described in Peterson et 
al. (1998). Additional improvements as suggested by Welsh 
(1999) are summarised in P04, which we use for our analysis. 

First an RSS procedure was applied by randomly se- 
lecting 22 observations from the hght curve. The flux uncer- 
tainties of the multiply selected observations were weighted 
according to Welsh (1999). This light curve was given as 
input to the FR procedure, where each measured value of 
the fluxes are modifled by adding the measured flux un- 
certainties multiplied with a random Gaussian value. The 
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Figure 4. The cross correlation function of the continuum at 5100 
A and the H/3 lightcurves, as obtained using the ICCF method, is 
plotted cLS ciS cl thick solid line. The auto correlation functions of 
the continuum at 5100 A and the Hp line are plotted as dashed 
and dotted lines respectively. The cross correlation function ob- 
tained using the DCF method is plotted as a thin solid line. 
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Figure 5. Histogram of the cross-correlation centroids obtained 
using the FR/RSS realisations. The solid and dotted lines repre- 
sent the cross-correlation centroid distributions (CCCDs) deter- 
mined using the ICCF and DCF respectively, for a bin size of 2 
days. 



Table 2. Estimates of centroid values obtained using the DCF and 
methods. 
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modified light curves were then cross-correlated and the cen- 
troids were determined as outlined above using CCF values 
above 50% of the peak value. This procedure was repeated 
for 4000 times, retaining only those CCFs whose maximum 
cross-correlation coefficient is large enough such that the 
correlation is significant with a confidence level of 95% or 
larger. A cross-correlation centroid distribution (CCCD) was 
built using the above centroids and is shown in Fig. [5] The 
average value of CCCD was taken to be Tcent- Since the 
CCCD is non-Gaussian (cf. Peterson et al. 1998), the upper 
and lower uncertainties in Tcent were determined such that 
15.87% of CCCD reahzations have r > Tcent + Ar„p and 
15.87% reahzations have t < Tcent — Atiow This error in 
Tcent corresponds to ±la errors for a Gaussian distribution. 

The centroid time lags obtained using different cross 
correlation methods for different bin sizes are given in Ta- 
ble [2] The variations due to different bin sizes are within 
the error bars. The DCF method gives a mean time lag 
of about 4 days, whereas the ICCF method gives a mean 
value of about 3 days. This difference is of the order of the 
estimated errors on the time lags. Considering these uncer- 
tainties, both the DCF and ICCF methods give time delays 
that are consistent with each other. This suggests that the 



^stimated time lag is not a spurious result for the sampling 
of the light curves presented here, and the results of the 
ICCF methods are reliable. To be conservative, for further 
calculations we use the time lag with the largest scatter cor- 
res ponding to a bin size of 3.5 days, obtained using the ICCF 
method. 

Based on this analysis the average observed frame time 
^Tag between the H/3 and the A5100 A continuum light-curves 
gjvas found to be 3.06lj g7 days. After correcting for the time 
J^flilation effects using the redshift of the source, we found a 
time lag of S.Ollj g4 days in the rest frame of the source. 

The wavelength coverage of our observations, also in- 
clude the Hq, line, and a time lag using the Hq line can also 
be estimated. However, by repeating the analysis procedure 
presented here, a reliable time lag using Ha line could not be 
found, because the correlation curves are too noisy. This may 
be due to the shorter duration of the observations and the 
relative fiux calibration procedure using the [O in] A5007 A 
line situated much farther apart in wavelength from H^ be- 
ing unreliable (Grier et al. 2008). Unfortunately, the nearby 
doublet [S II] AA 6716,6731 A narrow lines, which could be 
used for the relative calibration of H^, are too weak. Thus 
we do not estimate the time lag using the H^ line. 



3.3 Line width 

In order to relate the time lag to the mass of the black hole, 
an estimate of the line width, the line dispersion aiine, of the 
broad emission component of Ha, is required. Following P04, 
it is relatively straight forward and more practical to mea- 
sure oiine, the second moment of the profile, directly from 
the root mean square (rms) spectrum. Indeed, in the rms 
spectrum the constant components, or those that vary on 
timescales much longer than the duration of the observation 
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Figure 6. A zoomed up version of ttie mean spectrum (top panel) and corresponding rms spectrum (bottom panel). The flux is in the 
unit of 10^^'' erg cm~-^ A~^. The rms spectrum shows negligible contribution from the two [O III] narrow lines. 



vanish, thus largely obviating the problem of de-blending the 
lines. To obtain the rms spectrum, all the observed spectra 
were combined using the formula: 



5(A) = 



N 



(2) 



In Fig. [6] we show the root mean square (rms) spectrum. 
It can be seen that the two narrow [O in] lines have almost 
completely disappeared in the rms spectrum. 

The mean value of aune corrected for the instrumental 
response of the spectrograph and the associated uncertainty 
were obtained following the bootstrap method described in 
P04. From our observed 22 observed spectra, we randomly 
selected 22 spectra, irrespective of whether a particular spec- 
trum has already been selected or not. Since some of the 
spectra were selected multiple times, the mean value of the 
resultant number of spectra were smaller by 8. These ran- 
domly selected spectra were then used to construct an rms 
spectra from which aune was measured and corrected for 
the instrumental resolution of the spectrograph. This proce- 
dure was repeated 10000 times and the mean and standard 
deviation of these realisations are taken as aune and its un- 
certainty respectively. A distribution of cr;i„e values obtained 
from the bootstrap method is also shown in Fig[71 We thus 
estimate a line dispersion of = 1725 ± 105 km sec^^. 



3.4 Rblr and the mass of the black hole 

Using the rest frame time delay, the radius of the BLR is 
estimated to be Rblr = 2.53t'j';gg x 10"^ pc. In Fig. [8] we 
show the measurement of Rblr. and -^L^^^^^^ luminosity of 
the source presented here along with the most updated data 
set given by Bentz et al. (2009b) and the additional source, 
Mrk 290, given by Denney et al. (2010). The solid line is the 
relation obtained by Bentz et al. (2009b). From this figure it 
can be seen that our measurement lags are in agreement 
with the known Rblr — L relationship. 



The mass of the black hole was estimated using the 
formula in P04: 



Mbh = / 



Re 



G 



(3) 



where /S.V is the width of the line and G is the gravita- 
tional constant. The parameter / is a scaling factor, which 
takes into account the geometry and kinematics of the BLR. 
Onken et al. (2004) found an empirical value of / = 5.5, 
using a sample of AGNs having both reverberation based 
black hole masses and host galaxy bulge velocity dispersion 
(cr*) estimates. This value relies on the assumption that 
both AGNs and quiescent galaxies follow the same M-bh- 
(T* relationship (Ferrarese & Merritt 2000; Gebhardt et al. 
2000). For this particular scaling, the appropriate velocity 
width /SV is the line dispersion in the rms spectrum Gune 
(Bentz et al. 2008). Adopting the Onken et al. (2004) scal- 
ing factor and the (June measured from our observations, 
we estimated the mass of the SMBH in H 0507+164 to be 
Mbh = 9.62+° ?^ x 10® Mq. 



4 CONCLUSION 

We present for the first time monitoring observations of the 
X-ray source and Sy 1.5 galaxy H 0507-1-164, spanning a 
time period of about one month. We have obtained 22 nights 
of spectra during this period, with a mean sampling time 
of about 1.6 days. We measured an observed frame time 
lag of about 3.06lli'j''g7 days between the changes in the H/j 
emission line flux and the changes in the continuum flux at 
5100 A. After correcting for the redshift, we find a corre- 



sponding time lag of 3.01 



+0.42 



days in the rest frame of the 



source. From this measured time lag we deduced a size for 
the BLR of Rblr = 2.53t'i';55 10"^ parsec and estimated 
a black hole mass of 9.621^-73 x 10*^ Mq. Our estimate of 
RsLj? using the measured lag of is in agreement with the 
RsLii — AL^^^^^ relationship shown by Bentz et al. (2009b). 
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Figure 7. Histogram of tlie estimates of the crKne using the boot- 
strap method described in the text. 




Figure 8. The radius of the BLR vs tlic continuum luminosity at 
5100 A. The open circles axe from Bentz et al. (2009b). The open 
triangle is the source Mrk 290 reported by Denney et al. 2010. 
The solid line (with a slope of 0.519) is the fit obtained by Bentz 
et al. (2009b). The new measurement of this work is shown as a 
filled circle. 
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